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Abstract 

Background: Gene expression in a cell entails random reaction events occurring over disparate time scales. Thus, 
molecular noise that often results in phenotypic and population-dynamic consequences sets a fundamental limit 
to biochemical signaling. While there have been numerous studies correlating the architecture of cellular reaction 
networks with noise tolerance, only a limited effort has been made to understand the dynamic role of protein- 
protein interactions. 

Results: We have developed a fully stochastic model for the positive feedback control of a single gene, as well 
as a pair of genes (toggle switch), integrating quantitative results from previous in vivo and in vitro studies. 
In particular, we explicitly account for the fast binding-unbinding kinetics among proteins, RNA polymerases, 
and the promoter/operator sequences of DNA. We find that the overall noise-level is reduced and the frequency 
content of the noise is dramatically shifted to the physiologically irrelevant high-frequency regime in the presence 
of protein dimerization. This is independent of the choice of monomer or dimer as transcription factor and persists 
throughout the multiple model topologies considered. For the toggle switch, we additionally find that the presence 
of a protein dimer, either homodimer or heterodimer, may significantly reduce its random switching rate. Hence, 
the dimer promotes the robust function of bistable switches by preventing the uninduced (induced) state from 
randomly being induced (uninduced). 

Conclusions: The specific binding between regulatory proteins provides a buffer that may prevent the propagation 
of fluctuations in genetic activity. The capacity of the buffer is a non-monotonic function of association-dissociation 
rates. Since the protein oligomerization per se does not require extra protein components to be expressed, it 
provides a basis for the rapid control of intrinsic or extrinsic noise. The stabilization of regulatory circuits and 
epigenetic memory in general is of direct implications to organism fitness. Our results also suggest possible avenues 
for the design of synthetic gene circuits with tunable robustness for a wide range of engineering purposes. 



Background 

Recent experiments on isogenic populations of mi- 
crobes with single-cell resolution [?,?,?] have demon- 
strated that stochastic fluctuations, or noise, can 



override genetic and environmental determinism. In 
fact, the presence of noise may significantly affect the 
fitness of an organism [?] . The traditional approach 
for modeling the process of molecular synthesis and 
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degradation inside a cell is by deterministic rate 
equations, where the continuous change of arbitrar- 
ily small fractions of molecules is controlled instan- 
taneously and frequently represented through sig- 
moidal dose-response relations. However, the rate- 
equation approaches can not explain the observed 
phcnotypic variability in an isogenic population in 
stable environments. In particular, when molecules 
involved in feedback control exist in low copy num- 
bers, noise may give rise to significant cell-to-cell 
variation as many regulatory events are triggered by 
molecules with very low copy numbers < 100 [?]. A 
well known example is the regulation of inorganic 
trace elements [?], such as iron, copper, and zinc. 
While these trace elements are essential for the activ- 
ity of multiple enzymes, their presence may quickly 
turn cytotoxic unless their concentrations are care- 
fully controlled. 

Although the presence of phenotypic variation 
due to stochastic fluctuations need not be detrimen- 
tal for a population of cells [?] , elaborate regulatory 
mechanisms have evolved to attenuate noise [?] . Sev- 
eral systems-biology studies have recently focused 
on a select set gene-regulatory circuits, in particular 
those with feedback control. Feedback control cir- 
cuits have been identified as important for multiple 
species and proven responsible for noise reduction 
and increased functional stability in many house- 
keeping genes through negative autoregulation [?], 
long cascades of ultrasensitive signaling [?], bacte- 
rial chemotaxis [?] , and the circadian clock [?] . Ad- 
ditionally, recent studies on iron homeostasis [?,?] in 
E. coli highlight the noise-reducing capability medi- 
ated by small RNAs. 

Here, we study reversible protein-protein bind- 
ing as a novel source for genetic noise control. In 
particular, we have quantitatively analyzed the ef- 
fects of protein oligomerization on noise in positive 
autoregulatory circuits as well as a simple toggle- 
switch [?]. The all-or-none threshold behavior of 
positive-feedback circuits typically improves robust- 
ness against "leaky" switching. However, due to 
their functional purposes, gene circuits involved in 
developmental processes or stress responses that of- 
ten accompany genome- wide changes in gene expres- 
sion are intrinsically noisier than negative feedback 
circuits. 

It is frequently observed that transcription fac- 
tors exist in oligomeric form [?], and protein 
oligomerization is an important subset of protein- 
protein interactions, constituting a recurring theme 



in enzymatic proteins as well as regulatory proteins. 
Well studied examples include the A-phage repressor, 
ACI (dimer), the TrpR (dimer), LacR (tetramer), 
and Lrp (hcxadecamer or octamer). While many 
of the RNA-binding proteins dimerize exclusively 
in the cytosol, the LexA repressor [?], the leucine- 
zipper activator [?,?], and the Arc repressor [?] have 
been shown to form an oligomer cither in the cytosol 
("dimer path") or on the DNA by sequential bind- 
ing ("monomer path"). Previously, the efficacy of 
monomer and dimer transcription-regulation paths 
to reduce noise was separately studied for a negative- 
feedback autoregulatory circuit [?]. In contrast, we 
have focused on oligomerization in positive-feedback 
autoregulatory circuits, as well as genetic toggle 
switches based on the mutual repression of genes [?] . 
We find that cytosolic transcription-factor oligomer- 
ization acts as a significant buffer for abundance- 
fluctuations in the monomer, overall reducing noise 
in the circuit. Additionally, the noise-power spectral 
density is shifted from the low- to the high-frequency 
regime. In the toggle switch, cytosolic oligomeriza- 
tion may significantly stabilize the functional state 
of the circuit. This is especially evident for het- 
erodimerization. 

Yet another interesting case of ligand-binding- 
mediated receptor oligomerization has been re- 
ported [?,?], where the formation of various struc- 
tures of oligomers may act to buffer the intracellu- 
lar signaling against noise. Although our modeling 
and analysis is based on prokaryotic cells, we expect 
our main findings to be organism independent since 
protein oligomers, especially homodimers, is such a 
common occurrence across the species [?], with ho- 
modimers comprising 12.6% of the high-fidelity hu- 
man proteome [?,?]. 



Results and Discussion 

Dimerization breaks long-time noise correlations 
in autogenous circuit 

To evaluate the dynamic effects of protein-protein 
binding in positive-autoregulation gene circuits, we 
construct several alternative models of positive auto- 
genous circuits. Each model emphasizes a different 
combination of possible feedback mechanisms, and 
the network topologies considered can be grouped 
into the two classes of monomer-only (MO) and 
dimer-allowed (DA) circuits, according to the avail- 
ability of a protein-dimer state (color coding in 
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Fig. 1). Wc further group the DA circuits into 
three variations, DAI through DA3, depending on 
which form of the protein is the functional tran- 
scription factor (TF) and where the dimerization 
occurs. For DAI, we only allow the dimer to bind 
with the DNA-operator sequence (dimeric transcrip- 
tion factor, DTF), while for DA2 dimerization oc- 
curs through sequential binding of monomers on the 
DNA. In DA3, the protein-DNA binding kinetics is 
the same as in the MO circuit, hence monomeric 
transcription factor (MTF), with the addition of a 
cytosolic protein dimer state. While we will only 
present results for DAI in this paper, there is no 
significant difference for DA2 and DA3 [Additional 
file 1]. 

Note that the feedback loop is not explicit in 
Fig. 1 but implicitly included through the depen- 
dence of RNAp-promoter binding equilibrium on 
the binding status of the TF-operator pair. The 
sign (positive or negative) and strength of the feed- 
back control is determined by the relative magni- 
tude of the dissociation constants between RNAp 
and DNA which is either free or TF-bound. For 
instance, topology DAI has positive feedback con- 
trol if K 30 = k 30 /q 30 > K 32 = k 32 /q 32 , and K 30 
corresponds to the level of constitutive transcrip- 
tion (transcription initiation in the absence of bound 
transcription factor). For each topology, we study 
the dependence of noise characteristics on the ki- 
netic rates by varying the dimer lifetime, binding 
affinity, and the individual association/dissociation 
rates (see Table and Fig. 1). While we only discuss 
positive feedback control of the autogenous circuit in 
this paper, we have obtained corresponding results 
for negative feedback control [Additional file 1]. 

Fig. 2 shows a sample of ten representative time 
courses for the protein abundance. The effect of 
stochastic fluctuations is marked in the MO circuit. 
However, in all the DA circuits where the protein 
may form a cytosolic dimer we observe a significantly 
reduced level of noise in the monomer abundance. 
The suppression of fluctuations persists throughout 
the range of kinetic parameters that (so far) is known 
to be physiologically relevant (see Table 1). 

Calculating the steady-state distribution for the 
monomer and dimer abundances (Fig. 3) we ob- 
serve a clear trend that the monomer Fano factor 
(variance-to-mean ratio) is reduced as the binding 
equilibrium is shifted towards the dimer. This trend 
is conserved for all the investigated DA topologies 
(see Supplementary Information) . As long as dimer- 



ization is allowed in the cytosol, the fast-binding 
equilibrium absorbs long-time fluctuations stemming 
from bursty synthesis or decay of the monomer. 
When a random fluctuation brings about a sudden 
change in the monomer copy number, dimerization 
provides a buffering pool that absorbs the sudden 
change. Otherwise, random bursts in the monomer 
abundance will propagate to the transcriptional ac- 
tivity of the promoter, leading to erratic control of 
protein expression. It should be emphasized that 
this has nothing to do with the sign of regulation 
and is in agreement with the observations of Ref. [?] 
for negative autoregulation. Surprisingly, the mag- 
nitude of noise reduction in the positive autoregula- 
tory circuit is nearly the same as that for negative 
autoregulation which is typically considered a highly 
stable construct [Additional file 1]. 

A heuristic explanation can be found from Jaco- 
bian analysis of a deterministic dynamical system, 
which is justified for small perturbations around a 
steady state. When a random fluctuation shifts the 
monomer copy number away from its steady-state 
value, the decay toward the steady state can be de- 
scribed by the system Jacobian. The disparity in 
the magnitude of the (negative) eigenvalues of the 
Jacobian matrix for the MO versus the DA circuits 
signifies that the perturbed state is buffered by fast 
settlement of the monomer-dimer equilibrium. This 
buffering occurs before random fluctuation can ac- 
cumulate, possibly with catastrophic physiological 
effects, explaining the coarse long-time patterns ob- 
served in the MO model in contrast with the DA 
circuits (Fig. 2). 

Frequency-selective whitening of Brownian noise 

The dimerization process itself generates stochastic 
fluctuations on a short time scale. However, since 
this time scale is essentially separated from that of 
monomer synthesis and decay (orders of magnitude 
faster), dimerization effectively mitigates monomer- 
level fluctuations. The frequency content of the 
fluctuations is best studied by an analysis of the 
power spectral density (PSD), which is defined as 
the Fourier transform of the autocorrelation func- 
tion [?], originally introduced for signal processing. 
Fig. 4 shows the noise power spectra of DAI, and 
the distinction between the MO circuit and the DA 
topology is immediately evident. In particular, wc 
note the following two features, (i) A power-law de- 
cay with increasing frequency and (ii) a horizontal 
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plateau for the DA circuits. The power-law feature 
is explained by the "random walk" nature of protein 
synthesis and decay: The power-law exponent is ap- 
proximately 2, which is reminiscent of Brownian mo- 
tion (a Wiener process) in the limit of large molec- 
ular copy numbers. Compared to other commonly 
observed signals, such as white (uncorrelated) noise 
or l/f noise, protein synthesis/decay has a longer 
correlation time. If the autocorrelation function of a 
time course is characterized by a single exponential 
decay, as is the case for Brownian noise, the PSD is 
given by a Lorentzian profile, and thus, well approxi- 
mated by an inverse-square law in the low-frequency 
regime. We do not observe a saturation value for the 
MO circuit, and it is likely not in the frequency win- 
dow of physiological interest. This may especially be 
the case for circuits where the correlation times are 
long. 

The noise reduction is in the physiologically rele- 
vant low- frequency regime, and in Fig. 4 we have in- 
dicated the typical values for a cell cycle and mRNA 
lifetime. Although stochastic fluctuations impose a 
fundamental limit in cellular information processing, 
multiple noise sources may affect cellular physiology 
non-additively. For a living cell, fluctuations are es- 
pecially relevant when their correlation time is com- 
parable to, or longer than, the cell cycle. At the 
same time, short-time scale fluctuations (relative to 
the cell cycle) are more easily attenuated or do not 
propagate [?] . Additionally, the observed flat region 
in the PSD of the DA circuits implies that as far as 
mid-range frequency fluctuations are concerned, we 
can safely approximate them as a white noise. This 
insight may shed light on the reliability of approxi- 
mation schemes for effective stochastic dynamics in 
protein-only models. 

Increased lifetime of dimer plays an important role 

The virtue of the cytosolic dimer state is also directly 
related to the extended lifetime of proteins when in 
a complex. Except for the degradation tagging for 
active proteolysis, a much slower turnover of protein 
oligomers is the norm. This is partly explained by 
the common observation that monomers have largely 
unfolded structures, which are prone to be target of 
proteolysis [?]. It has also been pointed out that 
the prolonged lifetime of the oligomeric form is a 
critical factor for enhancing the feasible parameter 
ranges of gene circuits [?] . As seen from Fig. 3 (also 
Table 2), the fold change of the noise reduction, 



while still significant, is not as strong for the (hy- 
pothetical) case of dimer lifetime being the same as 
that of the monomer (72/71 = 1/2). However, the 
low-frequency power spectra still exhibit almost an 
order-of-magnitude smaller noise power than in the 
MO circuit with the same rate parameters (Fig. 4). 
Hence, the noise reduction capability holds good as 
long as the dimer lifetime is kept sufficiently long 
compared with the monomer-dimer transition. 

Effects of homo-dimerization in genetic toggle 
switch 

The exceptionally stable lysogeny of the phage A, for 
which the spontaneous loss rate is < 10~ 7 per cell 
per generation [?,?], has motivated the synthesis of 
a genetic toggle switch [?]. Toggle switch is con- 
structed from a pair of genes, which we will denote 
as gene A and B, that transcriptionally repress each 
other's expression. This mutual negative regulation 
can be considered an effective positive feedback loop 
and provides the basis for the multiple steady states. 
The existence of multistability, in turn, may be ex- 
ploited as a device for epigenetic memory or for de- 
cision making [?]. 

As the general attributes of positive feedback 
with cooperativity suggest, a genetic toggle switch 
responds to external cues in an ultrasensitive way: 
When the strength of a signal approaches a thresh- 
old value, the gene expression state can be flipped by 
a small change in the signal. For example, the con- 
centration of protein A (B) may rapidly switch from 
high to low and vice versa. However, previous stud- 
ies of a synthetic toggle switch have shown that the 
noise-induced state switching is a rare event [?,?,?]. 
In the ensuing analysis, we aim to delineate the ori- 
gin of this exceptional stability. 

In a simple model, the monomer-only (MO) tog- 
gle, regulatory proteins only exist in monomeric 
form. Although an external signal is not explic- 
itly included, random fluctuations in the abundance 
of the circuit's molecular components will occasion- 
ally flip the toggle-state for the two protein species. 
Drawing on the results from our analysis of pos- 
itive autoregulatory gene circuits, we hypothesize 
that dimerization in the regulatory proteins of the 
toggle switch will serve to stabilize its performance 
against noise. We allow the protein products of each 
gene to form a homodimer, being either AA or BB, 
which is similar to the cl-cro system in phage A [?] . 
The dissociation constant for the dimers is defined 
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as K\ — qi/ki, where k\ is the rate of two monomers 
forming a complex, and qi the rate of the complex 
breaking up into its two constituents. 

We evaluate the effect of the fast protein binding- 
unbinding dynamics on the toggle switch perfor- 
mance by using either (i) the monomers or (ii) the 
homodimers as the functional form of the repressor. 
Fig. 5 shows, for selected values of the dissociation 
constant Ki, representative time series of the pro- 
tein monomer (left) and dimer (right) abundances 
for the case of (a) monomeric or (b) dimeric tran- 
scription factors, respectively. A careful analysis of 
the phase space (in presence of noise) for our chosen 
set of parameters confirms that the studied toggle- 
switch systems are in the bistable region [?]. 

When monomer is the functional form of the re- 
pressor molecule (Fig. 5(a)) and K\ is large (limit 
of low dimer affinity), the protein populations are 
dominated by monomers. Hence, the circuit effec- 
tively behaves as an MO toggle. As K\ decreases, 
we see that the level of random switching is sup- 
pressed: Analogous to the autogenous circuit, the 
dimer pool stabilizes the protein monomer popula- 
tion. However, the noise suppression is not mono- 
tonic with increasing dimer binding affinity Indeed, 
for very large binding affinities (small K\), the num- 
ber of random switching events is increased since 
the monomer is only available in low copy num- 
bers. Consequently in this limit, it becomes more 
likely that a small fluctuation in the monomer abun- 
dance can cause a dramatic change in the overall 
gene expression profile. The noise-stabilizing effect 
of dimerization is also reflected in the corresponding 
PSDs [Additional file 1]. For instance, we observe a 
marked suppression of low-frequency fluctuations in 
the monomer abundance with increasing K\. 

In Fig. 5(b) we show corresponding sample time 
series for the case of a dimeric repressor, all other 
properties being the same as in (a). While the over- 
all trends are similar, we do note the following dif- 
ference. Contrary to the monomeric repressor case, 
there are very few toggle events in the strong binding 
limit: Since the signaling molecules (dimers) of the 
dominant gene (the "on" -gene) tend to exist in large 
copy numbers, a significant fluctuation is needed to 
flip the state of the toggle switch. In the case of 
monomeric repression, the signaling molecule exists 
in low abundance in this limit. Thus, the domi- 
nant protein species in the dimeric-repressor system 
is able to maintain much better control over the state 
of the toggle switch. 



In Fig. 6, we show the distribution (Na — Nb), 
the difference in molecule abundance for the two 
protein species in the case of monomeric (left) and 
dimeric (right) transcription factor. The asymme- 
try with respect to the zero axis is caused by our 
choice of initial conditions (protein species A in 
high concentration and species B in low concentra- 
tion), as well as the finite length of the time se- 
ries. For monomeric transcription, the presence of 
dimers with moderate binding affinity sharpens the 
monomer abundance distribution while accentuating 
its bimodal character. This is in agreement with the 
qualitative observation from Fig. 5 on switching sta- 
bility. For dimeric transcription, we clearly observe 
that the symmetry of the system is broken for small 
values of K\, indicating that the state of the toggle 
switch is extremely stable, and hence, likely deter- 
mined by the choice of the initial conditions. 

To systematically quantify our observations on 
the interplay between dimer-binding affinity and the 
functional stability of the toggle switch, we gener- 
ated long time series (« 3 • I0 7 sec) to measure 
the average spontaneous switching rate. In Fig. 7, 
we show the average toggle frequency relative to 
that of the MO toggle for the binding affinities 
ifi/nM = {2,20,100,1000}, and the average MO 
switching rate is 7.5 x 10 _6 /hour. As expected, we 
find that intermediate values of K\ are able to sta- 
bilize the toggle switch. Fig. 7 also highlights the 
increased stability of the toggle switch for a dimeric 
versus monomeric transcription factor, the dimeric 
switching rates always being lower and approaching 
zero for strong dimer binding. 

Heterodimerization in genetic toggle switch 

We have also considered the case of heterodimeriza- 
tion in the toggle switch, since the noise- and func- 
tional stabilization of the switch may be directly af- 
fected by the composition and source of the dimers. 
Note that, the gene-regulation activity is conferred 
by the two monomer proteins A and B and not the 
heterodimer AB. However, we find that the pres- 
ence of (inactive) heterodimers gives rise to very 
similar noise-stabilizing effects as that of homod- 
imers (Fig. 7). In fact, the existence of heterodimer 
state allows the dominant protein species to effec- 
tively suppress the (active) monomers of the minor- 
ity species. Thus the heterodimer circuit shows dra- 
matically enhanced functional stability as compared 
to the case of homodimeric repressors, not sharing 
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the discussed vulnerability of MO circuit to intrinsic 
noise. Although, to our knowledge, this is a purely 
hypothetical toggle-switch design, it provides a gen- 
eral strategy for noise control in synthetic gene cir- 
cuits, along with previously proposed approach of 
overlapping upstream regulatory domains [?]. 



Conclusions 

Cells have evolved distinct strategies to combat the 
fundamental limits imposed by intrinsic and envi- 
ronmental fluctuations. We investigated the role of 
protein oligomerization on noise originating from the 
random occurrence of reaction events and the dis- 
crete nature of molecules. Recent efforts to correlate 
network structure with functional aspects may pro- 
vide valuable insights into approaches for network- 
level noise control [?]. While negative feedback is 
one of the most abundantly observed patterns to 
achieve the goal of stability, it begs the question 
of how cells reliably change the expression of genes 
from one state to another. The ultrasensitive re- 
sponse circuit, exemplified by the ubiquitous signal 
transduction cascades in eukaryotic cells, has been 
proposed as an answer to this question [?,?]. 

In addition to the combinatorial expansion of 
functional specificity, we argue that the availabil- 
ity of oligomeric states contributes to the attenua- 
tion of stochastic fluctuations in protein abundance. 
In positive autoregulatory gene circuits, where the 
abundance of an expressed protein controls its own 
synthesis rate, dimerization provides a buffer serv- 
ing to mitigate random fluctuations associated with 
the bursty transcription-translation process. We 
find that short-time binding-unbinding dynamics re- 
duce the overall noise level by converting potentially 
pathological low-frequency noise to physiologically 
unimportant, and easily attenuated, high-frequency 
noise [?]. 

Noise-induced switching generally signals a de- 
fect in cellular information processing. Untimely 
exit from latency in the lambda-phage system di- 
rectly implies, as the immediate consequence to 
viruses, increased chance of being targeted by a host 
immune system. In the case of a bacterium, the ex- 
pression of a specific set of sugar uptake genes when 
the sugar is absent from the external medium is a 
considerable waste of cellular resources. For exam- 
ple, lac operon of E. coli can be considered to have 
the circuitry of mutual antagonism between the lad 



gene and lactose uptake-catabolic genes [?]. A dif- 
ference lies in the non-transcriptional deactivation of 
the allosteric transcription factor Lacl. LacY, lactose 
permease, indirectly regulates Lacl by increasing lac- 
tose uptake, which in turn catalytically deactivates 
Lacl. Likewise, many pili operons of Gram-negative 
bacteria are also known to utilize heritable expres- 
sion states, which are of crucial role in pathogene- 
sis [?,?]. 

We expect that the random flipping of gene ex- 
pression states in the examples of positive-feedback- 
based genetic switches may very well be closely cou- 
pled with the fitness of an organism. Phenomeno- 
logical models relating the fitness of an organism to 
random phenotypic switching in fluctuating environ- 
ments have provided important insights into the role 
of noise [?], but still many questions remain unan- 
swered. 

Applying these insights to the design of a syn- 
thetic gene switch demonstrates the potential use 
of affinity-manipulation for synthetic biology, where 
the construction of genetic circuits with tunable 
noise-resistance is of central importance. In par- 
ticular, our analysis highlights the potential util- 
ity of heterodimerization to stabilize ultrasensitive 
switches against random fluctuations. In practice, 
small ligand molecules may be employed to regu- 
late and tune the binding affinity of regulatory pro- 
teins, being either monomers or dimers. Our results 
further suggest that the structure of the protein- 
interaction network [?] may provide important in- 
sights on methods for genome-level noise control in 
synthetic and natural systems. 



Methods 

Model construction 

To evaluate the general role of protein oligomeriza- 
tion in a broad functional context, we studied the 
two most common motifs found in genetic regulatory 
circuits: positive autoregulation and the bistable 
switch. The reaction scheme studied is summarized 
in Fig. 1, where the binding/unbinding reactions be- 
tween RNAp and promoter or between TF and oper- 
ator are made explicit. Each distinct binding status 
of DNA is associated with a unique transcription 
initiation rate, and then the overall rate of mRNA 
synthesis is a weighted average of the initiation rates 
for distinct binding status, where the weights are 
given by the relative abundance of each configura- 
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tion at equilibrium, determined by the calculation 
of binding energy [?]. Note that, neither binding 
equilibrium nor empirical Hill-type cooperativity is 
assumed ad hoc. In particular, we split the lumped 
transcription process into two separate events, (i) 
isomerization of closed RNAp-promoter complex to 
its open form and (ii) transcription elongation fol- 
lowed by termination. This is to reflect the avail- 
ability of the free promoter while the transcription 
machinery proceeds along the coding sequence of a 
gene as soon as the promoter region is cleared of the 
RNAp holoenzyme. Otherwise, the promoter would 
be inaccessible during a whole transcription event, 
altering the random mRNA synthesis dynamics. 



To realize the genetic toggle switch in a stochas- 
tic setting, we keep track of the microscopic origin 
of cooperativity that gives rise to bistability. Among 
various strategies, we employ multiple operator sites 
which have the same binding affinity with the re- 
pressor. The resultant circuitry is, in essence, two 
autogenous circuits, A and B, which are connected 
through the active form of their expressed proteins 
(the active form being either monomer or dimer). 
The connection is implemented by allowing the ac- 
tive form of proteins A (B) to bind the operator 
sites of gene B (A). In order to make the interac- 
tion between the two genes repressive, unlike positive 
autogenous circuits, K31 and K 32 in Fig. 1 arc now 
greater than K 30 , making the protein transcriptional 
repressors. For reasons of analytical simplicity, we 
have studied the symmetric toggle switch, where 
the reaction descriptions of each component follow 
those of the autogenous circuit. Again, the quan- 
titative characteristics of macromolecular binding- 
unbinding are chosen based on the phage lambda- 
E. coli system. The only exception is related to 
the multiple operator sites, where the second re- 
pressor binds an operator site with higher binding 
affinity when the first site is already occupied by 
the repressor protein [?]. We introduce three differ- 
ent dimerization schemes. Three different dimeriza- 
tion scheme have been introduced: (i) homodimer- 
ization with monomeric repressor, (ii) homodimer- 
ization with dimeric repressor, and (iii) heterodimcr- 
ization with monomeric repressor. By solving for the 
stationary states of the deterministic rate equations, 
we could identify the bistability region in parameter 
space to which all the model systems under consid- 
eration belong. 



Stochastic simulation 

While the deterministic rate equation approach 
or Langevin dynamics explicitly gives the time- 
evolution of molecular concentration in the form 
of ordinary differential equations, chemical master 
equation (CME) describes the evolution of a molec- 
ular number state as a continuous-time jump Markov 
process. To generate the statistically correct trajec- 
tories dictated by CME, we used the Gillespie di- 
rect [?] and Next Reaction (Gibson-Bruck) [?] al- 
gorithms, both based on the exact chemical mas- 
ter equation. The Dizzy package [?] were used as 
the core engine of the simulations. To ensure that 
calculations were undertaken in a steady state, we 
solved the deterministic set of equations for steady 
state using every combination of parameters inves- 
tigated. We employed these deterministic steady- 
state solutions as initial conditions for the stochastic 
simulations. For each model system, we generated 
10 5 ensemble runs with identical initial conditions 
and used the instantaneous protein copy number at 
a fixed time point t = 5000 sec. To achieve high- 
quality power spectra in the low- and high-frequency 
limits, we ran time courses (~ 10 5 sec) with higher 
sampling frequency (20 measure points per sec). 

To calculate the average switching rate, we gen- 
erated time series of minimum length 3 • 10 7 sec (ap- 
proximately corresponding to 1 year). We identify 
a state change in the toggle switch by monitoring 
the ratio of the monomer and dimer abundance for 
the two protein species. In order to avoid counting 
short-time fluctuations that do not correspond to a 
prolonged change of the toggle state, we a applied 
sliding- window average to the time series, using a 
window size of 1000 sec. 
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Tables 



Table 1 - Probability rates for positive autogenous circuit. 
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Reaction 


Value (s- 1 ) 


Rcf. 


protein dimerization 


ki 

9i 


Pi 
P 2 


+ 


Pi- 
Pi 4 


>P 2 
Pi 


0.001-0.1 
0.1-1 


[?,?] 




feo 


P 2 


+ 


D00 


-> D20 


0.012 






920 


D20 - 


- P 2 


+ D00 


0.9 




TF-operator int 


k 2 i 


Pi 


+ 


D00 


-> D10 


0.038 




921 


D10 - 


- Pi 


+ D00 


0.3 


[?,?,?] 




^22 


Pi 


+ 


D10 


-> D20 


0.011 






922 


D20 - 


-» Pi 


+ D10 


0.9 







&30 


R + 


D00 


D01 


0.038 






930 


D01 


-» R 


+ D00 


0.3 




RNAp-promoter int 


&31 


R + 


D10 - 


-+ Dll 


0.038*, 0.38* 




931 


Dll 


-> R 


+ D10 


0.3*, 0.03* 


[?,?,?] 




k 3 2 


R + 


D20 - 


D21 


0.38** 






932 


D21 


-» R 


+ D20 


0.03** 




Isomerization 


V 


Dxl 


-► C 


+ DxO 


0.0078 


[?] 




a 


C -> 


M + 


R 


0.03 






P 


M - 


> Pi 4 


- M 


0.044 




tsx-tsl elongation & decay 


7o 


M -> 




0.0039 


[?,?] 




7i 


Pi - 


+ 




7xl0~ 4 






72 


P 2 - 


+ 




0.7-3. 5xl0" 4 





Kinetic rates for the positive autogenous circuit. Experimentally available rates are all taken from lambda 
phage- E. coli complex. The values with superscript correspond to the circuit topologies DAI (*), DA2 (f), 
and DA3 (f) in Fig. 1. 



Table 2 - Relative Fano factors of protein abundance distributions 



Ki (nM) 


72 = 7i/10 


72 = 7i/2 


monomer dimer 


monomer dimer 


1 


0.127 0.809 


0.132 0.679 


20 


0.209 0.936 


0.230 0.716 


500 


0.866 0.478 


0.826 0.426 



The Fano factor of protein abundance distribution for the autogenous circuits (topology DAI), relative to 
that of the monomer-only (MO) circuit, 8.729. 
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Figures 



Figure 1. Schematic of model autoregulation gene circuit. 

The DNA binding status is indicated by Dxy, where x corresponds to the operator region (empty=0, 
monomer=l, dimer=2), and y to the promoter region (empty=0, RNA polymerase bound=l). C repre- 
sents the open complex of DNA-RNAp holoenzyme with the promoter sequence just cleared of RNAp and 
is subject to transcription elongation. Finally, M, PI and P2 correspond to mRNA, protein monomer, and 
dimer, respectively. The network topologies can be grouped into two classes, monomer-only (MO) or dimer- 
allowcd (DA) circuits. We have studied DAI (red lines), which only allows the dimer to bind with the 
DNA-operator sequence, DA2 (green) with sequential binding of monomers on the DNA, and DA3 (blue), 
which shares protein-DNA binding kinetics with MO while allowing dimerization in the cytosol. Note that 
for topology DA2, we have chosen K 3i = K 30 (see text for details) We have assumed cells to be in the 
exponential growth phase and the number of RNAp (R) constant. 



Figure 2. Ten independent time courses of the abundance of protein monomers in the (positive) 
autoregulatory circuit. 

The availability of a cytosolic dimer state (red, using circuit topology DAI) significantly reduces the copy- 
number fluctuations of the monomer compared to the monomer-only (MO) circuit (blue). All corresponding 
MO and DAI parameters have the same values. In the ensuing simulations initial conditions are chosen to 
be the steady state solution of the corresponding deterministic rate equation so that the transient behavior 
should be minimized. 



Figure 3. Stationary state distribution of monomer (black) and dimer (orange) protein abun- 
dance in the positive autogenous circuits. 

The left (right) column corresponds to a ratio of the dimer and monomer decay rates of 72/71 = 1/10 
(72/71 = 1/2). The molecular copy numbers are collected at a fixed time interval (5 • 10 3 sec) after the 
steady state has been reached. Here K\ = q\jk\ is the dissociation constant of the protein dimer. As 
the binding equilibrium is shifted towards the dimer state (decreasing K\), the noise level is monotonically 
reduced (see Table 2). Note that the prolonged protein lifetime due to the complex formation (left column) 
affects the noise level. 



Figure 4. Power spectral density (PSD) of fluctuations in protein abundance. 

The PSD of the MO circuit clearly displays a power-law behavior. All other model systems with an available 
cytosolic protein dimer state (DAI shown here) develop a plateau in the mid-frequency region regardless of 
the model details (see Supplementary Information). As the dimer binding affinity increases, the noise level 
is further reduced. We have included the MO result in the dimer panel (right) for reference. Datasets with 
solid (empty) symbols correspond to 72/71 = 1/10 (72/71 = 1/2). 



Figure 5. Sample time series of monomer and dimer copy numbers in genetic toggle switch. 

(a) MTF circuit, where monomer is the functional form of the repressor, (b) DTF circuit, where dimer is the 
functional form of the repressor. The left (right) column shows the number of the two monomer molecules 
A and B (dimers AA and BB) , and the initial state is always with species A (red) in high abundance. Note 
that the switching frequency depends on the binding affinity of protein dimer. 
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Figure 6. Distribution of monomer abundance differences between protein species A and B. 

The asymmetry with respect to the zero axis is due to the choice of initial state (species A high) and the 
finite time span of simulations. 



Figure 7. Random switching rates of genetic toggle switches. 

Ordinate is the ratio of the random switching rates of various toggle switches to that of the monomer-only 
(MO) circuit, 7.5 x 10~ 6 /hour. MTF, monomeric transcription factor; DTF, dimeric transcription factor; 
Het-MTF, monomeric transcription factor with deactivated heterodimer state. 
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Additional Files 



Additional file 1. Supplementary results for the positive autoregulatory circuits with various 
topology. 

Protein abundance distribution and power spectral density of autogenous DA2 and DA3 circuits are pre- 
sented. 
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long dimer half-life 



short dimer half- 




Number of molecules 



Monomer 
Dimer 



Dimeric repres 



1 - K = 2 nM 

1 




A B 



A B 



1.2 



1.0 - 



0.8 - 



■2 0.6 - 



0.4 - 



0.2 - 



0.0 - 



— i 

K.= 2nM 



Monomer, MTF 

Dimer, MTF 

- Monomer, DTF 

- Dimer, DTF 
Monomer, Het-MTF 




T 



K = 20 nM 

1 



T 



K = 100 nM 

1 



